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ABSTRACT 

Research on wind and temperature profiles over land have 
established a number of fundamental relationships that are tested 
here over an oceanic location. Most attempts at calculating the 
roughness parameter over oceanic surfaces have been based on 
surface-layer theories, and often employ an iterative procedure in 
approaching the sea Surface. This paper employs a similarity 
relationship between the wind and temperature profiles, using data 
at ship "Victor" appropriate to forced convective theory, as 
defined by Priestley. 

The theory of forced convection, as applied here, makes it 
necessary to compute the Richardson number and the Monin-Obukhov 
scale length as preliminaries to the computation of roughness 
parameter. The roughness parameter and friction velocity are 
obtained from the theory using an iterative procedure. The wind 
speed at 20 meters is correlated against scale length and roughness 
parameter. The latter makes very little contribution to the variance 
of the windspeed while the scale length contributes in a manner 


reasonably expected from theory. 
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List of Symbols 


Symbol Definition 
un, Wind speed at n meters height 
Z Height in centimeters above the 


water surface 


k Von Karman: Constant 

uN Friction velocity 

L Monin-Obukhov scale length 

T Scaling temperature 

il Eddy fhgat transport in the. 
vertical 

Rj Richardson number 

g Acceleration of gravity (980 cms) 

Cy Specific heat of air 

P Density of air 

Zo Roughness parameter 

4 Potential temperature 

S Non-dimensional wind shear after 


Monin & Obukhov 
Ky Eddy diffusivity for heat conduction 


Ku Eddy diffusivity for momentum flux 


THitwordTreerion 5 


K6-20 Ratio of wind at 600 cms to that 
at 2,000 cms. 
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Symbol Definition 
Ky 4 Ratio of z #4 )y 4 to the wind 
: difference at 20 meters minus that 
atten centimeters 
T Surface stress 


A more exact form of K) 4 


Vg/fZ, Surface Rossby number 





1. Introduction 

Studies of the roughness parameter, Zo, have been numerous 
over land and a number of fundamental relationships using various 
meteorological parameters have been determined. Far less 
research has been accomplished in this field over the oceans. 

This is readily understandable when one considers the expense and 
difficulty in erecting and maintaining an instrumented " platform 
at sea. With this in mind, one could turn to what should be a 
prime source of data--the ocean station weather ships. 

Certain problems become evident very quickly when 
considering station Victor. Four different ships occupy this 
station and only two have the same anemometer height. This, at 
first, seemed to present an insurmountable problem in using the 
data but the theoretical error involved in using a median height of 
twenty meters instead of actual heights is less than two and one 
half percent. When considering that winds are recorded only to 
the nearest knot and temperature to the nearest degree, the error 
introduced in the wind profiles by considering all heights equal to 
twenty meters is insignificant. 

Monin and Obukhov [3] , in their universal similarity theory 
introduced three scaling parameters that are considered essentially 


invariant with height for layers near the ground. These three 
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parameters are scaling velocity, u*, scaling length, L, and 
scaling temperature, ee and have been well described in a 
recent text (Lumley and Panofsky) [J on atmospheric turbulence. 


ive u* and T* are defined as follows: 


Ls 
P= Mee (1) 


P=4 5 
=~ FG P (2) 


va —— (3) 


It should be noted that a selection of cases has been limited 
to cases for which -Ri@.032 holds. According to Priestley's drifi- 
rion(47J our restriction limits the study to the cases of "forced 


convection". Here Ri refers to the turbulent parameter 


Ri = gi pez 


7\2 (4) 


and is to be distinguished from the flux Richardson number Rf, 


defined as 


Re — Ke RF; (5) 





Under normal circumstances, it is considerably easier to deal 

with Ri than with Rf, since the factor Kj; / Kyy is itself an 

increasing function of Ri. 

2. Important properties of the forced-convective surface layer 
According to the similarity theory for the near-neutral, or 

forced convection case (as defined in section 1), the wind 


profile can be written 


¥ Iz Z 
U= $0,342) (6) 
In a like manner the temperature profile can be written, after 


_Lumley and Panafsky 
6-Q=T (nE+ €F) (7) 


In both (6) and (7), Lumley-and:Panof sky [1] Suggest the value 


/ 
eC = 4.5. 
By dividing equation (7) by equation (6), and taking the 


derivative of both numerator and denominator with respect to 


height, one has a relationship governing the constants T* and 


u*/k: 





_ — “su/az (8) 
) 





Relationship (8) means that the change in the temperature 
over a given height is proportional to the change in the wind over 
the same height. Even more important, this ratio is constant 
with height since tT, uu”, and k are all considered constant in the 
surface layer. 

The ratio of these constants rf (u* /k) may be determined 
from the air temperature at psychrometer height (six meters) and 
the temperature at Zp , assumed to be the same as the sea 
surface temperature, so that equation (8) may be written, in 


finite-difference form centered at three meters as 
Fe 
( ) — Us-Uo 


Here u, is the wind speed at Z, and is zero. By equation (8) this 
ratio is constant at all levels in the surface layer, including that 
at 1.41 meters. 

From the definition of L, it follows that L satisfies the 


identity 


Z= ft kis (10) 


However, Priestley /f4Jhas introduced a surface layer scale 


lenght L' given by 


~——ae or (1) 





so that (10) becomes 


f=RS (12) 


a formula lacking the troublesome ratio Kyz/Kyy which has been 


found to depend upon Ri. But from the definition of T*, we have 


Ger = Gp 
Ke Gp laZ) = Ku gleaZ (13) 


From (13), together with the similarity profiles (6) and (7), which 


Or 


apply in the case of forced convection, it follows that Ry=Kyy 
and therefore L=L'. This is a fundamental property characteristic 
of similarity profiles in forced convection. 

3. Derivation of a specific value for (Ri)) 4] 


From equations (4) and (9), the Richardson number is 


expressible in the form 


ies i ro (Z3H72),, (14) 


where every parameter in (14) is now known except for the 


quantities ug and (73 4° 





The surface~layer wind profile is assumed to be valid to 
heights of twenty meters. The windspeed Ug may then be related 


to Ugg by a new parameter through 


Us = Uae Ke-2e (15) 


where Ke-20 is formed from the ratio of Ug to UgQ: This now 
leaves the quantity (Z oy y i as the sole unknown in 
equation (14). One may approximate this quantity by forming the 
ratio of two wind differences centered at 1.41 meters. In the 
following equations, subscripts denote the height of a variable 
in meters although substitution in centimeters is generally 


required. For example we have 


(Z 52 jh : RA rs Hg (16) 


and 


U ~ UL, 
Peat 2% 
U Cos ie, -— 


AZ OO 1 on 
Therefore, it follows that 


ea = Ur20 Ky we) 


Now the last two unknowns in equation (14) are found in terms of 


(17) 





known or measured quantities. Following Martin PJ 


K.4 = = (i (19) 
ee #8) 





and, similarly, 


, oe o + es 
Ke-zo = | oO ET. 


with all values in centimeters. Ke_99 may be approximated 


(20) 


wRpNs 


satisfactorily for obtaining Ri by omitting In Zo. 
The right sides d (15) and (16) are now known as a function of L 


so that (14) may be solved for Ri as 


Lie O — (Q-@) 
HI » K (Le) Koo o(L) O, wh 6 (21) 


where K, 4, and Kg-20 are given as functions of Lin (19), (20). 
In the forced convection regime, Lumley and Panotsky show 


that the Richardson number may be written as 


A i = 
i -_ | + eZ ) fos = Km (22) 
4 L, 
Eliminating Ri between (21) and (22) leads to a quadratic in L, 


the only meaningful root of which is given by 


2 7 ve 


A = 1.49, 200)-.7hn'Y5-9000 Ri (h, 2 00+4n 2000) 
B 


) I 


24 Ri (dn 200\(Ln 2.000) (24) 


C = Ie Rin 200Ln 2000\81x10") 
~ 400 4 Ri (Ln 2 00M L. 2000)] 


where 





Ri may now be obtained by solving equation (22) with Z=141.1 
cm and the value of L obtained from (23). 

The quadratic equation defined by (23) and (24) had as 
input data only windspeeds at 20 meters, and temperature - 
differences measured from the surface to 6m. Only one data-set 
was selected at ship Victor per day, and these were selected with 
a view toward insuring cases of forced convection: for example, 
a vertical temperature difference less than or equal to 6T., and 
windspeeds as large as possible under these circumstances were 
selected. A table of Ri values computed as a function of both UxQ 


and 6g - © is included in the Appendix. 





4. Iterative procedure of solution of roughness parameter 
Following MartinL 2 7], use is made of equation (6) to obtain 


the quantity Ugg ~ U9] in the form 
“es x 
Uzo~ Us = (2.302 6)(2.30103)+ 2, (ApO0 Lr") (25) 
But with Uy9 as follows 


Uno = 2.302 6 (Ht) dope 22 + (H' Me) 26) 


one obtains the result, Up, ;, given in (27) 


Uo = 005 Upo + 48 (2.2988 XI-BS4H5 Loo Zo) 27 


By using the Monin-Obukhov normalized wind shear 
,_ Rk 7 OU 
> Ut & OL. 


and (12), it follows that 


f= (Bea) (OP), MAB), 


Because Z, is small compared to 20 meters, the quantity 
( dU/a Ln Z } 4 will be rewritten involving a larger span:, as 
follows 


eke = Ci 


(iol = 





where ©) 4 will be defined immediately below. 


AU vs - = 
AT Z ), = eS = Cry (Uzo~ Up.i) 


From (6), we then.have 


4 644 


Ci - Mpabh /(Uzo -Uai) = aoq } + (29) 


Combining (28) and (29) leads to 


f= BE), Ca ( Uru) 


where C, 4 is defined in (29). 


A third equation used in the iteration-process is 


Or Ul (2.3026) - bra Zo \\+ fe | (31) 


but this has already been built into equation (25). The latter may 
be used for more accurate determination of Z, once an earlier 
iterate, or satisfactory guess, has been established. Equation 
(31) assumes the use Of equation (6) to an elevation of 10 cm above 
mean sea level. 

To adopt an efficient iterative method which does not lay 
too much emphasis on the last 10 cm, we must give proper 
weight to the layer 0-20 m. Here, we use equations (27) and (30) 


with ug _; eliminated. This leads to the result 
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U* — 1.4 (LE!) (0.495) 


~ Cin (GB). (22488-2211 eyez) (32) 


and, of course, we also have 


U20 = (Le 200 + bn 2 ae 1000 ) (33) 


The iteration process may be summarized as follows: 

1. Assume fey = 0.1 and obtain a corresponding estimate 
of (u* /ky from (32). 
(1) 


2. With this first estimate of (u*/k)‘-’, use equation (33), 


with ugg and L known, to determine Za 
3. Repeat the procedure using Z._{2) in (32) to get a second 


estimate of (u* /k) (2) ; 


4. Repeat the iteration until a final . residue is obtained: 


J Zo" ~ 25 


When this has been accomplished bye is taken to be 2 , and (i/e\”’ 
is the appropriate value from either (32) or (33). 

One safeguard must be followed in the use of (32): the 
denominator of the right side must not be allowed to become zero 


or negative. However, the use of (33) helps to restrict (u* /k) to 


reasonable, non-negative values. 


ll 





The results, giving Z,, u'/k as well as ugg, Ri’, Land Ri 
are listed in Table 1 (Appendix). In addition, some statistical 
relationships between the parameters are also considered in the 
following section. The raw data considered was drawn from 30 
randomly selected dates. The value of Lis computed from 
equation (23), while Ri, 4] is obtained from (22). 

The computations were made using data cards obtained from 
the National Weather Records Center, Asheville, North Carolina, 
for ocean station "Victor". All cards selected were for 0300Z, the 
observation nearest noon local time, from the months of November 
and December of 1960 and of March 1961. The number of 
intervals between successive days was considered a random 
function of the date. Wind values were tabulated in this paper in 
centimeters per second, and temperatures in degreesCelsius. Sea 
conditions during the period were not considered, with the 
exception that no major storms were in the immediate area. 

5. Resultssandsconclushemssion s 

The values of u99, Land ay for the sample of 30 cases were 

subjected to a multiple regression analysis leading toa result of 


form 


Un, = -SOHEZ-O0053IL+75797 — (34) 


variable 3 variable 1 variable 2 
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Of the two independent variables tested in (34), Z, plays an 
insignificant part in that its partial regression coefficient was 
-0.00969. When the mean value of Ze = 0.1066 cm was inserted 


into (34), the simple regression equation 
(35) 

Uso = -0.0053IL +75 4.76 
resulted. The variable L had a partial correlation coefficient of 
-0.56153, so that its effect upon the reduction of the variance 
of Ugg was sufficiently strong that resulting F-value due to the 
regression involving L was 6.576. This F-value, with 2 and 27 
degrees of freedom, is Significant at a level in excess of the 99 
percent confidence level. This result is not surprising from the 
manner in which 1/L enters U20, Eq. 33. A brief summary of some of 


the other statistics obtained are listed: 


Ub, 218707 ens =, = 0.10662 cm L=-43,773.bem 


Cy 


37847ems Ge, =0.10258em T= 40,682.bem 
6.2270.0096! gy 20.5615 3 
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The most surprising fact regarding the statistics, and 
equation (34) is the complete unimportance of the variable Zo in 
describing u99. At first glance, equation (33) would suggest a 
much stronger dependence, since -1nZ, is one of the terms of the 
right side, and this should indicate that ugg decreases with 
increasing Z,5. This is a result presented by many writers, who 
emphasize the tendency for u to decrease with increasing surface 
Rossby number Vg/fZ,5. However, from equation (33) it can be 


seen that 


(36) 





where (u*/k), is the neutral value of (u”/k). That part of (u*/k) 


which depends upon stability is primarily incorporated in the ratio 


(|+ ete até) /(j+ |e aun 


so that the neutral part of (u* /k) is essentially 


(4) 7 Cu 


The final point to be emphasized here then is as follows: 
Within equation (36) the term -inZ> tends to reduce ug9. However, 


carefully measured results of Hay (1955) [5] indicate that 
2 
Ue = 5.14 Zo 
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so that u may be expected to increase in the boundary layer as Z, 
increases. Obviously these two effects Oppose one another, with 
a wind-increasing tendency from (u*/k) and a wind-decreasing 
tendency arising from -1nZ, within the bracket of (36). This 
oddity resulting from opposite effects produced by an increase in 
Zo appears to explain the lack of correlation of u79 with Zo. 

No statistical analysis involving u*/k with Zo has been made 
but their values are listed in the Appendix, and may serve as the 


basis of future research on the subject. 


RS) 
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APPENDIX I 


Table of selected micrometeorological parameters in the surface layer. 





a At(°c) rit (10°9) L Ri u%/k Z. 
1081.5 - 3.88 mis -25197.329 -.0057 1a, Sale .098703 
412.0 -4.44 -9.02 -4900.552 -.0327 56.065 .250295 
309.0 -1.66 -6.00 -6399.332 -.0243 oe 457 - 34632 
1030.0 -3.88 -1.26 -23040.770 -.0063 107.509 . 060164 
1339.0 -2.22 -0.43 ~64212.942 -.0022 171.079 . 709676 
1081.5 = 27 7 -0.82 -34474.195 -.0041 111.440 01238 
1442.0 -2.77 -0.46 -59822.933 -.0024 146.968 . 145697 
BLS. 0 -1.11 -1.44 -20410.493 -.0071 54.025 . 126881 
1442.0 -6.11 -1.01 -28242.996 -.0050 149.440 .06137 
1133.0 -2.77 -0.75 - 37639.751 - .0038 116.493 .042618 
618.0 = 222 -1.98 -15405.752 -.0095 65.760 . 095000 
1133.0 -1.11 -0.29 -92272.089 -.0015 114.868 .111753 
2060.0 222 -0.18  -150915.860. -.0009 208.067. .158484 
S755 -4.44 -1..98 -15426.860 -.0095 93.154 . 089827 
Hee 5 - .55 -0. 32 -85614.909 -.0017 78. 378 .098703 
618.0 -1.11 -0.99 -28764.428 -.0049 64.009 . 035126 
ee - 3.88 -2.55 -12423.038 -.0119 77.801 . 208904 
875.5 -1.66 -0.74 -37875.237 =. 0038 55.825 .174716 
1133.0 -2.79 -0.74 -37980.775 -.0038 116.468 .098703 
1030.0 -2.22 -0.71 - 39240. 355 -.0036 105.800 .098703 
721.0 -3.33 -2.18 -14163.219 -.0103 98.529 . 738132 
1133.0 -0.56 -0.15 ~-%84229.295 -.0008 104. 304 . 062660 
33.0 -3.89 -1.02 -29043.613 -.0051 117.4444  .098703 
1287.5 -2.78 -0.56 -49075.979 -.0029 147.914 B57 
1081.5 -1.11 -0. 32 -85030. 331 -.0017 109.737 .018613 
1236.0 -1.11 -0.25 +110239.246 -.0013 15a hs .057789 
1545.0 eel -0.86 - 32759.215 -,0044 148.858 .087884 
772.5 -1.11 -0.63 -44052.566 - .0032 79.153 .098703 
618.0 -1.67 -1.48 -19944. 649 -.0073 37.485 . 023398 
721.0 -3.89 =2.54 -12462.568 -.0118 65.775 . 103646 
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